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Abstract. A phase transition for bosonic atoms in a two-dimensional anisotropic optical lattice is consid- 
ered. If the tunnelling rates in two directions are different, the system can undergo a transition between 
a two-dimensional superfluid and a one-dimensional Mott insulating array of strongly coupled tubes. The 
connection to other lattice models is exploited in order to better understand the phase transition. Critical 
properties are obtained using quantum Monte Carlo calculations. These critical properties are related to 
correlation properties of the bosons and a criterion for commensurate filling is established. 

PACS. 05.30.Jp , 03.75.Lm, 67.90.-Hz 

1 Introduction esting many-body systems may be probed using a variety 

of optical techniques. 

Atoms in optical lattices offer new versatile ways of study- The optical potentials can be controlled with great 
ing many-body phenomena [HH]- Interference patterns of range and precision, giving access to any kind of Bravais 
laser light provide periodic potentials for atoms, which are lattice as well as superlattices and quasiperiodic potentials 
cooled to the nanokelvin range using laser and evaporative [3] . The geometrical setup and choice of polarisation de- 
cooling. In this way, ground-state properties of these inter- termine the lattice geometry, and the laser irradiance de- 
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termines the amplitude of the potential. If the laser light is 
detuned to the red of the atomic transition, the potential 
will typically form an array of wells in which the atoms 
may be trapped if the potential is deep enough. The height 
and width of the potential barriers separating the wells de- 
termine the rate of tunnelling between them. In this way, 
the dimensionality of the sample can be controlled: mak- 
ing the potential barriers high enough in one Cartesian 
direction will effectively inhibit all tunnelling and result 
in a stack of independent, two-dimensional (2D) lattices. 
Increasing the potential along a second direction yields a 
2D array of ID tubes of atoms. 

Recently, such dimensional crossovers in optical lat- 
tices and their relation to other theoretical models have 
attracted considerable attention. References [HEUS] stud- 
ied a 2D array of one-dimensional tubes of atoms, known 
as Tomonaga-Luttinger liquids (TLL) [7J. In a random- 
phase approximation, phase boundaries, coherence prop- 
erties and excitation spectra were expressed in terms of 
the parameters that characterise the TLLs. In Ref. [S], 
we studied the corresponding situation in a 2D lattice ge- 
ometry and investigated the relation between this model 
and a classical XY model [Oj . It was determined that the 
transition between independent TLL tubes and a 2D su- 
perfluid is of Berezinskii-Kosterlitz-Thouless (BKT) type 
[TUIITTIIT^ if the density is held fixed, and as an example, 
the critical point was determined for a specific value of 
tunnelling matrix element and filling. 

In this paper, we further investigate the transitions in 
the 2D case. Using quantum Monte Carlo calculations, we 
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investigate the link between coherence properties and the 
critical point for the phase transition. We confirm that the 
transition is determined by a TLL theory and that it de- 
pends crucially on the commensurability of the atoms in 
the lattice. In Sec. [2l we lay out the theory. In Sec. [3l we 
show how the present Hamiltonian is equivalent to other 
lattice models, which clarifies the nature of the phase tran- 
sition. Section [3] explains the numerical method and the 
finite-size scaling performed in order to locate the phase 
transition, and presents the results for the phase transi- 
tion at unit filling. In Sec. [5] we show how it is linked to 
the coherence properties. In Sec. [SI it is illustrated how 
changing one of the parameters of the theory allows for 
a crossing a two subsequent phase transitions. In Sec. [71 
we investigate commensurability and density dependence, 
and finally in Sec. [8] we summarise and conclude. 

2 Phase transitions in a 2D Hubbard model 

Consider a one-component gas of bosonic atoms in a 3D 
lattice potential created by three standing waves at right 
angles, 

V{r) = Vox cos^ ki^x + Voy cos^ ki,y + Vqz cos^ ki^z, (1) 

where ki^ — 2t:/X is the angular wave number of the hght, 
and the potential heights Vox, Voy and Vqz are determined 
by the laser irradiance and detuning. In the tight-binding 
approximation, the many-body Hamiltonian is expanded 
in a basis of Wannier functions w„(r — r^), each of which is 
localised in a well centred at position r^, where i labels the 
wells and n is a band index. If the wells are tight enough. 
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and the atoms are cold enough, only the lowest-energy 
Wannicr function within each well, wolr — r^), needs to 
be taken into account; this is the lowest Bloch band. The 
gas is then described by the bosonic single-band Hubbard 
model fT3l. 



<ij> i i 

(2) 

where i is a placeholder for the three indices ix,iy, iz enu- 
merating the lattice sites in the Cartesian directions, and 
the sum subscripted <ij> runs over pairs of neighbouring 
sites. The parameter U is the on-site interaction strength, 
tij is the tunnelling matrix element for the barriers be- 
tween sites i and j, and /i is the chemical potential, ad- 
justed in the calculations to give the desired density of 
atoms. The on-site interaction strength is defined as 



U ^ 



J d'^r|wo(r — r. 



(3) 



where a is the s-wave scattering length of the atoms and m 
is the atomic mass. Since this optical lattice is created by 
three independent pairs of standing waves, the tunnelling 
matrix element tij in one particular term takes on one of 
three values. If the wells i and j are neighbours along the 
X direction, its value is 



2m 



y(r) 



woix + d,y,z), 
(4) 

where d = A/2 is the lattice spacing. The tunnelling ma- 
trix elements ty and tz are defined similarly. 

It turns out that the tunnelling matrix elements de- 
pend exponentially on the potential depth 13J. As a re- 
sult, one can easily make tx, ty and tz differ by orders of 



2D SF 




Fig. 1. Predicted phases of the anisotropic Hubbard model 
in two dimensions. "2D MI" and "2D SF" stand for the two- 
dimensional Mott insulating phase and superfluid phase, re- 
spectively. The phase "ID MI" is predicted to exist only when 
one of the side lengths of the system is finite. In this figure, the 
system extends infinitely along the y direction and is 4 sites 
wide in the x direction. The slight decrease of the correspond- 
ing phase boundary is not visible in the figure. 

magnitude by a judicious choice of the laser irradiances. 
Thus, Vqz can be chosen large enough that tz can be ne- 
glected and the sample is effectively 2D. In the following, 
we confine the discussion to the 2D Bose-Hubbard model 
(henceforth referred to as the Hubbard model for brevity). 
The phase diagram is sketched in Fig. [1] and we now de- 
scribe the general features. When both t^ and ty are small 
compared with U, and in addition the number of bosons is 
commensurate with the number of lattice sites, the sam- 
ple is Mott insulating (MI), with exponentially decaying 
phase correlations and suppressed fluctuations in the num- 
ber of particles per site. In the opposite limit, when t^ 
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and ty are of the same order as U or larger, the ground array of finite tubes lining up along the y direction and 
state is the 2D superfluid (2D SF) state, characterised by extending along the x direction. 

long-range phase coherence in 2D and a fluctuating nuni- Because of the finite excitation energy within the tubes, 
ber of particles per site. In one dimension with one parti- there can now exist a ID MI phase if ty is sufficiently 
cle per site, there is a quantum phase transition between small. By the same argument, the tunnelling along the z 
these phases at t/U — 0.3 14J. In higher dimensions, one direction can be neglected, as we assumed, 
may apply a mean-field approximation [9 , dictating that In a finite system, there is no true phase transition. In 
the transition occurs when the sum of all matrix elements practice, this is not a problem since the phase transition 
equals the ID value, i.e., is replaced by a crossover. In theory, it is always possi- 

ble to let the system extend indefinitely along the weakly 

= 0.3 (5) 

coupled y direction, while keeping it finite along the x 

• 11 1 • 1 r.T^ direction. There is then a true phase transition between 
m two dimensions. The exact, numerically obtained 2D re- 

, _p. Tiir 1- nii insulating and superfluid behaviour along the y direction, 

suit [15j deviates slightly from this mean-field result, i'ur- 

^ , , , . , , We call this the ID MI - 2D SF transition, 

thermore, quantum Monte Carlo calculations performed 

In Refs. [Hinil^ , the ID tubes were described using 



on finite systems tend to underestimate the critical value 

_ , , . . TLL theory. This theory describes a ID many-body sys- 

for the phase transition. 

tem, which is characterised by two parameters indepen- 

If one of the tunnelling matrix elements, say t^, is rel- 

dently of statistics and the detailed properties of the con- 

atively large and the other, ty, is small, there may in fact 

stituent particles, namely the sound velocity Vg and the 

under certain circumstances exist a phase where there is 

TLL parameter K, defined as 

superflow along one dimension but not along the other. 

Such a state will be called the ID MI phase. It turns out x ~ — (Q) 

Vs' 

that in an infinite 2D system, such a state is absent: if 

Here, vp — Ti'Kp/ [md), where p is the filling factor, i.e., the 

there is superflow in the x direction, then any nonzero 

number of bosons per site, and d is the lattice constant. 

tunnelling along the y direction will put the sample into 

The TLL parameter K determines the behaviour of the 

the 2D SF state [T^ (as a corollary, any finite tunnelling 

particle-particle correlations r(i,j), which obey the power 

along the z direction will result in a 3D SF state). How- 
law 

ever, as was first noted by Ho et al. and Gangardt et al. 

mnn,7n , . ^ihj) = {a\aj) oc |r, - rj|-i/(2if)^ (7) 

[Hlollo] (although applied to the 3D case), the situation 

changes when the sample is finite along the strongly cou- For the discussion of the correlations along different di- 
pled x direction. The system can then be thought of as an rections, we introduce the notation Fx and Kx for the 
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correlations within the strongly coupled tubes, as well as the whole 2D system is superfluid, it may not be valid. It 
Fy and Ky for describing correlations along the weakly was found in Refs. [4l[5] that 
coupled y direction. 

Ej^tyN^-"-^, (10) 

Since a larger implies greater coherence along the 
X direction, we expect to increase with tx/U . Further- and 

C U 

more, it is a known exact result that Ec — — , (11) 



2 < Kx < oo (8) where the exponent 



for a system of lattice bosons with only on-site interac- q,^ = . (12) 

tions [7]. If one further decreases the tunnelling beyond 

Thus, it is predicted that the behaviour of the particle- 

the point where — 2, one enters the MI phase. How- 
particle correlation function within the tubes determines 

ever, the quantitative relation between and is not 

the location of the critical point for the decoupling of the 

known in closed form. 

tubes. Here, is the number of sites that the lattice 
extends in the x direction, and Cq is a constant that can in 

3 Equivalence to other lattice models 

general not be obtained in closed form. The number-phase 

In order to obtain an expression for the transition point, Hamiltonian undergoes a SF-MI quantum phase transition 

the TLL tubes were in Refs. [1[S] treated as structureless ^t the critical point where 

sites by integrating out all degrees of freedom except a E-j 

— ^ — constant, (13) 

number operator Nj and a phase operator 0^ for each 

tube (the analysis was done with a 3D system in mind, where the right-hand side is a universal numerical constant 
but the results apply to the 2D case as well). In that way, to be determined. Hence, for the critical value tyc of ty, 
a number-phase model 1171 was obtained with a governing f 
Hamiltonian 

as long as the analysis of Ref. [HIS] provides a valid model 

H = -Ej^ cos{q^,-<P,) + Ec^iN,-Nof, (9) 

<^ijy i for the anisotropic optical lattice. The same functional 

where A^o is the equilibrium number of particles per tube, form for the critical point was obtained in Ref. 6 by 
and <ij> denotes neighbouring tubes in the y direction, means of a random-phase approximation. In Eq. (|14p . the 
Ej is usually referred to as the tunnelling energy and Eq filling p is defined as the mean number of atoms per lat- 
as the charging energy. The mapping is expected to hold tice site, and is thus dimensionless. Due to the known con- 
in the ID MI phase, but above the phase transition, where straints on K^, the power on lies between -1.75 and -2, 
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where the former holds close to the transition to a 2D MI. where Co is the same constant as in Eq. (fTTjl . and (3 = 

We define the effective coupling constant 1/(^3^), where T is the temperature of the Hubbard model. 

~ ty T2-a i-a /1c;^ This mapping shows that if the XY model has a phase 

= JJ^x (15) 

transition at an inverse temperature /?xy, then the Hub- 
Stated in terms of this quantity, the critical value of the 

bard model has a quantum phase transition when ty reaches 

tunnelling ty can be written as 

a critical value, consistent with Eq. (jl6p . The critical point 

jj ^yc^x P 1 V-^"7 will depend on the tunnelling within the tubes, ia;, through 

where and iy^ are constants to be determined. t^^c Luttinger parameter and also through the un- 

The nature of the quantum phase transition in the known constant of proportionahty Co in Eq. 
number-phase model can be understood by making use of The 2D XY model exhibits a BKT transition [TOllISl 

the general result that a quantum phase transition in D UHlin]- This kind of transition possesses several charac- 

dimensions is in the same universality class - i.e., it has teristics, among them a universal jump in the superfluid 

the same critical properties - as a classical phase transi- density at the transition point. According to the series of 

tion in a certain corresponding model in 13 + 1 dimensions mappings performed here, from a 2D anisotropic Hubbard 

[9]. In ID, with the number of bosons kept fixed, the tran- model via a ID number-phase model to a 2D XY model, 

sition studied here is in the same universality class as the the anisotropic Hubbard model exhibits a BKT transi- 

classical XY model in 2 dimensions at a finite tempera- tion. The findings of Ref. confirm these expectations, 

ture TxY = 1/(^b,9xy)- The mapping is accomplished by and the present paper further expands on the subject, 
identifying [TS] 

\ Ec' 4 Finite size scaling 

Lx,XY ~ EjEc, 

T T /I '7^ In the numerical calculations, the 2D Hubbard model was 

^a,XY = ^y, (17) 

simulated using the stochastic series expansion method 
P^23j . The chemical potential /i was tuned to obtain 
the desired mean number of particles. By selecting those 
Monte Carlo steps that correspond to a fixed number of 
particles, it was made sure that the calculations were per- 



where quantities with the subscript XY refer to the 2D XY 
model and quantities without that subscript refer to the 
number-phase model. In terms of the underlying anisotropic 
Hubbard model, we can thus write 



XY 



CqU J ' formed at a given density, which is important for the char- 

1/2 

Lx,XY = fi {CoUtyL^ ""p^ , acteristics of the phase transition. We chose f3 = lOOOf/"^ 

Ly^XY = Ly, (18) in order to make sure that ground-state properties were 
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calculated, and the number of states per site was chosen 
to 6 in order to ensure convergence. 

Systems of x Ly sites were simulated, where the 
side lengths and Ly were varied in order to assess the 
predicted dependence of the critical point on the length 
discussed in Sec. [5J The boundary conditions were chosen 
to be periodic, which is necessary for calculating the su- 
perfluid density as will be described. The most important 
calculated quantities are the superfluid density ps and the 
particle-particle correlations and Fy. The superfluid 
density is numerically computed via the winding number 



Wy as 



(19) 



where Wy is the net number of times that a particle line 
crosses the periodic boundary in the y direction in the 
simulations [10 . The superfluid density as a function of 
the ratio of tunnelling matrix element and on-site interac- 
tion energy, ty/U^ is plotted in Fig. [2] for a few examples 
of parameter values. It is clear that the curves make a 
sharp drop at a transition point when ty/U becomes small 
enough. The transition point seems to depend on both 
and p, as was predicted in Sec. [31 and we now turn to the 
problem of calculating this dependence. 

In order to locate the critical point, finite-size scaling 
must be performed [TTlll2l[Tn] . However, the analysis of 
the present problem presents several difhculties compared 
with a classical XY model. At the BKT transition, the 
superfluid density assumes a value known as the univer- 
sal jump, proportional to the critical temperature. This 
is routinely used in 2D XY- model simulations, but since 




Fig. 2. A few arbitrarily selected curves showing the superfluid 
density as a function of the ratio of tunnelling ty to on-site 
interaction U . The different curves correspond to different data 
sets obtained at different fixed side lengths Lx and Ly, and 
different filling factors p. 

in the present case the two quantities are only known 
to within an unknown constant, we cannot make use of 
this relation. Moreover, the mapping between the Hub- 
bard model and the 2D XY model is only valid at and 
below the phase transition. In the 2D SF phase, the co- 
herence between tubes may be comparable to that within 
the tubes and the mapping is not valid. Finally, a quan- 
tum Monte Carlo calculation of the 2D Hubbard model 
is very time consuming and in the parameter regime of 
interest, where the parameters tx and ty differ by orders 
of magnitude, it is hard to obtain data with high accu- 
racy. This is aggravated by the fact that the side length 
Ly varies between 4 and 32 in the simulations. 

As seen in Sec. [31 the anisotropic Hubbard model is 
predicted to have the same properties close to the transi- 
tion as the 2D XY model has close to the BKT transition. 
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if the 2D XY model has a superfluid density given by 



Ps.XY 



T 1 — CKa: 



(20) 



The BKT transition occurs when ty assumes a critical 
value; at this point, all the curves Ps,XY{ty, L^', Ly) com- 
puted for different parameter values should ideally coin- 
cide (taking finite-size effects into account). This means 
that has to be optimised so that all the curves for 
different parameter values coincide as closely as possible. 
This is accomplished by first considering sets of data series 
with a given Ly and different L^. ff the analysis in Sec. 
[3] is correct, then all systems of size Lx x Ly can be ap- 
proximated as ID chains of length Ly, and therefore the 
results for Ps,XY as a function of ty should coincide be- 
tween data series with similar Ly and different Lx, if only 
the parameter ax is chosen correctly. As noted in Sec. |3l 
the correspondence is only expected to hold in the f D MI 
phase, so that for ty > tyc, one cannot require the data to 
coincide. Thus, one needs to consider the variance among 
the curves below a supposed critical point, and choose the 
value of ax that minimises the variance. This is to be done 
for each Ly separately, and then the results for different 
Ly can be compared. We thus do a spline interpolation of 
the points over the relevant range of ty and compute, for 
a given Ly and ax, 



Y'Av{ps^XY){iy,ax,Ly) = 

^ ( Ps,XY{iy;Lx; Ly) - —— ^ 



Ps.xvity; Lx'jLy) 



0.01 
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Fig. 3. Variance in the superfluid density among data sets as 
a function of the unknown exponent = l/{2Kx)- 

and furthermore the summed variance over the whole range 



of ty is 



Tj-. . ^ fdiyV£iT{ps,xY){iy,ax,Ly) 
Va,T[ps,xY)[ax,Ly) J— J- . (22) 

In Eqs. (I2m22p . Nl^ is defined as the number of different 
values of Lx used. Since the curves are supposed to coin- 
cide only below the critical point, the integration limits for 
ty are chosen as < ty < 0.3, anticipating the result that 
the critical point is close to tyc = 0.3. The result does not 
depend strongly on the chosen integration limits. For each 
value of Ly, we obtain an optimal value for ax- Averaging 
over different Ly, the result is 



resulting in 



ax = 0.15, 



Kx = 3.4. 



(23) 



(24) 



(21) 



When data for different system sizes Ly are compared, 
they are expected to coincide at the critical point, but not 
below or above. In addition, the coincidence of the curves 
is exact only in the limit Ly ^ oo, but it is known how 
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Fig. 4. Superfluid density for the dual XY model, scaled using 
the best fit parameter Kx ~ 3.4, for a filling p = 1. The plot 
contains all the data sets used in the scaling. The values of Ly 
are as indicated in the legend. 

to make the lowest-order correction for finite Ly. At the 
BKT transition point in the 2D XY model, the superfluid 
density depends asymptotically on the size of the finite 
sample as 

Ps,XY(iy) 



Ps,xy(oo) 



1 



(25) 



2\n{Ly)+C 

where Ps,XY(iy) is the superfluid density computed us- 
ing a side length Ly. This result is known as Weber- 
Minnhagen scaling [TT] , and in Ref. [in| it was found that 
the procedure applies to non-quadratic systems as long 
as the winding number is computed along the shorter di- 
mension of the sample. In our effective XY model, L^^xy 
is proportional to the inverse temperature of the Hubbard 
model and it is thus much larger than Ly. The constant 
C was in Ref. [T2] found to be equal to 1.8. 

Figure S] collects the simulated data for all different 
values of L^ and Ly. Ideally, all the data is expected to 
coincide at the critical point tyc. The variance among these 



Table 1. Parameters determining the 2D SF-ID MI transi- 
tion. The method to calculate the error bars are described in 
the text. The data for tx/U = 0.3 are lifted from Ref. [§]. 



tx/U 


P 


tyc 


Otx 


Kx 
(transition) 


Kx 

(correlation) 


0.3 


1.0 


0.3 


0.25 


2.0 




0.5 


1.0 


0.33 


0.15 


3.4 


3.0 


0.5 


5/4 


0.32 


0.15 


3.4 


2.7 


0.5 


19/16 








2.1 



curves is now recorded as a function of ty and the smallest 



variance is at 



tyQ — 0.33. 



(26) 



Table [T] summarises the computed parameters for a few 
different values of tunnelling and flUing p. The con- 
stant tyc was defined in Eq. (jl6p , the exponent and the 
result for are as obtained above, and the value for 
obtained from correlations will be discussed in Sec. [5] 

In order to check the above results, we apply a dif- 
ferent scaling procedure, by bunching together results for 
the same L^ but different Ly. If the superfluid density, 
corrected as in Eq. ([25|l . is computed for a range of Ly 
values and a fixed L^, the dependence on L^ cancels out 
and the curves should coincide at the critical point. The 
dependence of the critical point on L^ can then be calcu- 
lated. Figure [5] shows how this method works. The point 
tyc/U, at which the variance of the Weber-Minnhagen 
scaled superfluid density across different Ly is a mini- 
mum, is recorded for each fixed value of L^. (The result 
for Lx — S is indicated with an arrow in the upper panel 
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Fig. 6. Particle-particle correlation function in the x direction. 
The correlation function is shown as a function of coordinate 
for filling p = 1, t:c/U = 0.5, = 32, Ly = 8, ty/U = 
2.32 X 10"" (*), and ty/U = 2.32 x 10"^ (o). 



itative support for the prediction of Refs. [UllllS] that the 
critical couphng decreases as a power-law function of L^;, 
with a power slightly below 2. 



Fig. 5. Upper: An example of data collapse for the super- 
fluid density scaled according to Weber-Minnhagen scaling, at 
side length Lx ~ 8. Lower: Measured critical tunnelling as a 
function of side length Lx; the dependence is expected to be 
a power law. The linear fit yields ax — 0.17, giving Kx ~ 2.9. 
The circle indicates the previously obtained critical point for 
the Mott transition in the ID Hubbard model, Lx ~ 1- 

of Fig. [5l) Then tyc/U as a function of is fitted to a 
power-law dependence, as illustrated in the lower panel of 
Fig. [5] The best linear fit to the log-log-curve is given by 
tc ~ 0.31L^^-^^. In the figure, we have also inserted the 
previously obtained result ty^. = 0.30 for the case = 1, 
which is just the ID Hubbard model [M]. Fig. [5] is a qual- 



5 Phase transition and correlations 

One important prediction made in Refs. [HISIIB] is that the 
dependence on the transition point tyc on tube length 
is linked to the behaviour of the particle-particle correla- 
tions in an isolated tube. In order to test this, we compute 
the correlation function Fxi^ix), where is the number of 
lattice sites separating two points in the x direction. An 
example of a computed correlation function is shown in 
Fig. [51 The correlation function is fitted to a power law 
according to Eq. ([7]). It is seen in Fig. [6] that the corre- 
lations in the strongly coupled x direction depend on the 
tunnelling in the weakly coupled y direction (just as the 
opposite relation holds). Since the predictions of Refs. [H 
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Fig. 7. Superfluid density in the x and y directions, respec- 
tively, for a Bose-Hubbard model of size 12x12 lattice sites, 
fixed tunnelling ty/U — 0.003, and a filling of one atom per 
site. The scale for ps,y has been expanded by a factor 200. The 
points labelled A, B, and C are discussed in the text. 

[51[5] build on the correlation properties of an isolated tube, 
we should use the results obtained for the smallest values 
of ty, in the ID MI phase, where the tubes are decou- 
pled. We find the value — 3.0 for — 0.5U and filling 
p = 1. This is consistent with the value = 3.4 found 
from finite-size scaling in Sec. |31 

6 Dependence on in-tube tunneling 

As a way to visualise the three predicted phases, we show 
as an example in Fig. [7] the result of a calculation where 
the tunnelling in the strongly coupled x direction, t^, has 
been changed while ty is kept constant at ty — 0.003C/. 
This corresponds to moving along a horizontal line in the 
lower part of the phase diagram in Fig. [TJ (Observe, how- 
ever, that Fig.[l]corresponds to the case — 4, while here 
Lx — 12, and therefore the position of the phase bound- 
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ary is shifted.) The parameters are chosen such that, an- 
ticipating the results in Sec. [H the system should pass 
from the 2D MI phase, via the ID MI phase, into the 2D 
SF phase. The calculation is done for a finite lattice with 
12x12 sites, and the superfluid densities corresponding to 
motion in the x and y directions, respectively, are found. 
It is seen that the superfluid density corresponding to the 
X direction begins to increase flrst, at the point labelled 

A, and the increase of ps.y commences at the later point 

B, This is precisely what is expected for a system that 
crosses the two transition lines. However, one should note 
that the increase of ps^y seems to saturate at point C, and 
at the same point, the slope of the curve for ps^x is also 
seen to slightly decrease. 

In fact, when one tunes t^ such that the 2D MI, ID 
MI, and 2D SF phases are visited in turn, then both phase 
transitions belong to the BKT universality class. Further- 
more, experience shows that it is the onset of the drop 
from a finite value, not the onset of a rise from a value close 
to zero, that should be identified with the BKT transition 
point. In Fig. [7] it thus seems that it is the point labelled 

C, rather than the points labelled A and B, that indicates 
the true transition, and one can conclude that the transi- 
tions associated with the strong and weak couplings occur 
(within calculated error bars) at the same point. This is 
also what the theory for the infinite system dictates: there 
is no ID MI phase if both and Ly are taken to infinity. 
Here, however, we are concerned with finite systems, so we 
should not take the limit of infinite L^- The figure shows 
clearly that in the flnite system, there is a region in which 
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the superflow along the x direction is non-negligible but 
that along the y direction is very small. This is the ID 
MI phase, but it cannot be found by applying finite-size 
scaling for the x direction. 

7 Phase transition at non-integer density 

We return to Eq. (fTH]) . obtained from Refs. [4 l [5 l l6] , which 
predicts that the transition point tyc/U has a power-law 
dependence on both the length and filling p of the 
tubes. However, the ID MI phase only exists if there is 
an integer number of bosons in each tube, i.e., if the fill- 
ing is commensurate with respect to the number of tubes. 
The number of particles per site, on the other hand, does 
not have to be an integer. To check this, and thus estab- 
lish that we are indeed seeing a Mott transition along one 
direction, we study the cases p — 5/4 and p — 19/16, re- 
spectively. In the first case, the tube length, L^, is chosen 
as a multiple of 4 in order to ensure commensurability, 
and in the second case, we choose values of that are 
not divisible by 16, in order to avoid commensurability. 
The results are shown in Fig. [51 

The finite-size scaling was performed as described in 
Sec.m In the commensurate case, p ~ 5/4, a phase transi- 
tion is found. For the incommensurate case, p — 19/16, the 
data may be collapsed with the best-fit result K,^ = 3.6, 
and iyc — 0.21. However, it is seen in Fig.[5]that the super- 
fluid density for this incommensurate density does not go 
steeply to zero when ty is decreased below the calculated 
critical point, even for the largest system size Ly — ^2. In- 
stead, the data is clearly consistent with the curves meet- 
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Fig. 8. Superfluid density in the XY model, ps,XY, for a num- 
ber of data sets. Data is taken at the commensurate density 
p = 1 (o); commensurate tube filling with p = 5/4 (v); and 
the incommensurate filling p = 19/16 (*). Here, tx/U = 0.5, 
the value of Lx for the difi'erent curves varies between 2 and 
32, and Ly between 4 and 32. The quantities on the axes are 
the best fits to finite-size scaling for the coupling and super- 
fluid density, respectively. The inset shows the same data in a 
log- log plot. 

ing the x axis at the origin, unlike the commensurate cases 
p — I and p ~ 5/4. We conclude that the data corrobo- 
rates the conclusion that the filling per tube needs to be 
integer in order for the ID MI phase to exist. 

8 Conclusions 

In this paper, we have studied bosons trapped in two- 
dimensional optical lattices by Monte-Carlo calculations, 
with the objective of characterising phase transitions and 
their dependence on dimensionality and lack of isotropy. 
Apart from the expected phases, where one have super- 
fluidity or Mott insulation along both directions, we show 
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